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ABSTRACT 


This  report  offers  an  introduction  to  the  transfer  matrix 
method  of  analyzing  the  dynamic  behavior  of  common  engineering 
structures,  followed  by  an  explanation  of  the  application  of 
the  transfer  matrix  method  to  an  array  of  aircraft  panels 
which  are  continuous  over  supporting  stringers.  The  skin- 
stringer  problem,  important  to  the  prediction  of  fatigue  fail¬ 
ures,  is  discussed  for  rather  general  conditions.  The  rect¬ 
angular  panels  may  vary  in  thickness  and  length  while  the 
stringers  may  vary  in  cross-sectional  shape  and  size.  The 
panels  may  be  flat  or  curved,  and  the  curved  panels  may  vary 
in  radius  of  curvature. 
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Air  Force  Materials  Laboratory,  Wright-Patterson  Air  Force  Base. 
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I*  INTRODUCTION 


The  following  presentation  is  an  attempt  to  clearly  outline  the 
transfer  matrix  method  of  analyzing  the  dynamic  behavior  of  an  elastic 
system  and,  at  the  same  time,  to  explain  a  recent  extension  of  this 
theory  to  aircraft  type  panel- stringer  construction^!  ,  For  additional 
introductory  information  the  reader  is  referred  to  Reference  [  2]  , 

The  technique  of  transfer  matrices  is  related  to  the  methods  of 
dynamic  analysis  developed  by  Holzer  t  ^  J  and  Myklestad!4!  *  Like 
these  previous  methods,  the  transfer  matrix  method  is  a  calculation 
of  the  deflections  and  internal  forces  at  successive  values  of  a  single 
space  coordinate  (stations)  by  utilizing  a  knowledge  of  the  system 
intertia,  damping,  and  stiffness  properties  between  stations.  Thus, 
the  similarity  of  these  procedures  extends  to  the  computation  of  sys¬ 
tem  natural  frequencies  by  iterative  procedures,  to  the  determination 
of  normal  modes,  and  to  the  calculation  of  deflection  and  force  and 
moment  type  responses  in  the  case  of  forced  vibrations.  In  other  words, 
the  transfer  matrix  method  accomplishes  the  same  types  of  objectives 
in  much  the  same  manner  as  the  Holzer  and  Myklestad  methods.  The 
difference  between  these  older  techniques  and  the  transfer  matrix 
technique  lies  in  the  advantage  of  conciseness  that  results  from  the 
use  of  matrix  algebra  by  the  latter  method.  This  advantage  has  made 
practicable  the  analysis  in  detail  of  complex  structures  such  as  single 
rows  of  curved  panels  supported  at  varying  intervale  by  not  necessarily 
symmetrical  stringers*  The  use  of  matrix  algebra  does  not  compro¬ 
mise  the  original  advantages  of  the  Holzer -Myklestad  style  of  analysis. 
For  example,  the  transfer  matrix  method  allows  the  introduction  of 
appropriate  and  separate  damping  descriptions  at  each  component  of 
the  structure.  Considering  the  advances  that  are  currently  being 
made  in  damping  technology!^!  ,  this  is  no  small  advantage-  On  the 
other  hand,  the  transfer  matrix  method  leads  to  certain  difficulties  in 
numerical  computation.  These  problems  and  their  remedies  that  have 
proven  to  be  effective  will  be  briefly  discussed. 

To  explain  the  principles  of  the  transfer  matrix  method  a  series 
of  examples  of  increasing  complexity  will  be  employed.  These  examples 
also  illustrate  the  scope  of  this  procedure  which  is  presently  limited  to 
a  structure  undergoing  sinusoidal  motion  (or  being  in  static  equilibrium-- 
a  special  case),  and  to  a  structure  being  one  dimensional  in  space  or 
whose  mathematical  description  is  reducible  to  that  where  the  unknown 
amplitude  factor  is  a  function  of  just  one  variable*  An  example  of  the 
latter  alternative  is  a  thin  plate  of  rectangular  shape  which  is  simply 
supported  at  two  opposite  edges.  In  the  case  of  a  single  frequency 
vibration,  we  could  write  for  the  vertical  deflection 


1 


w(x,y,t)  -  e^11  Yn  (y)  sin 

where  y  is  the  space  coordinate  parallel  to  the  simply  supported  edges, 
and  Yn(y)  describes  the  variation  of  w  in  the  y  direction.  In  other 
words,  we  have  reduced  a  plate  problem  to  a  one  dimensional  problem 
susceptible  of  solution  by  transfer  matrix  methods  by  (justifiably) 
assuming  the  form  of  the  deflection  with  respect  to  one  of  the  space 
coordinates  * 

The  examples  in  the  order  of  their  discussion  are  (1)  an  undamped 
beam  with  discrete  masses,  (2)  an  undamped  and  damped  beam  with  distrib¬ 
uted  mass  properties,  (3)  the  calculation  of  frequency  response  functions 
for  a  beam,  (4)  a  flat  panel-stringer  row,  and  (5)  a  curved  panel-stringer 
row-  The  first  three  examples  are  for  explanatory  purposes-  The  last  two 
are  summaries  of  recent  extensions  of  transfer  matrix  methods. 
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II •  the  myklestad  beam  problem 


A  convenient  place  to  begin  an  explanation  of  the  principles  of 
the  transfer  matrix  method  of  analysis  is  with  an  undamped  straight 
beam  constructed  of  discrete  masses  connected  by  massless  beam 
segments — the  Myklestad  "lumped  mass"  beam  model.  Figure  1  is 
a  drawing  of  a  general  Myklestad  model  beam,  where  y,  and  J1  are 
the  discrete  mass  and  rotary  inertia  at  station  j,  respectively.  We 
could  begin  our  analysis  by  examining  either  of  the  two  essential 
features  of  this  beam,  the  elastic  supports  and  the  concentrated  masses 
located  at  the  stations,  or  the  beam  segments  between  stations.  Let 
us  begin  by  taking  a  detailed  look  at  typical  massless  beam  segment  (j). 
Figure  2,  where  for  generality  we  will  include  all  possible  combinations 
of  stiffness  and  inertia  at  the  stations  that  mark  the  segment  endpoints. 

If  in  an  actual  application  any  of  these  properties  are  absent,  their 
descriptive  constant  is,  of  course,  zero.  The  combination  of  a  linear  and 
a  torsional  spring  could,  for  example,  be  supplied  by  another  beam  of 
small  mass  running  perpendicular  to  the  beam  under  consideration  and 
attached  to  it  at  the  given  station. 

To  locate  a  position  along  this  massless  beam  segment  we  will 
make  use  of  the  local  coordinate  Xj .  If  Wj (xj ,  t)  is  the  vertical 
deflection  of  the  segment,  positive  down,  the  equation  of  motion  is 
easily  found  from  the  basic  beam  equation  to  be 


(Ei)  ^Vj/ax.4  =  o  (i) 

where  for  the  sake  of  simplicity  we  have  chosen  the  beam  segment  to 
have  uniform  stiffness  properties.  (If  we  could  not  reasonably 
describe  the  stiffness  properties  between  our  desired  mass  stations 
as  being  uniform,  we  could  of  course  approximate  the  non-uniform 
stiffness  by  piecewise  uniform  segments  separated  by  stations  with¬ 
out  inertia  properties,  or  even  if  deemed  necessary  at  the  cost  of 
complicating  the  equation  of  motion,  we  could  describe  the  non- 
uniform  stiffness  as  an  analytic  function  of  xj.)  If  we  are  concerned 
with  a  single  frequency  vibration  of  the  beam,  we  know  that  every 
point  of  the  beam  will  move  in  a  sinusoidal  fashion;  i.e.,  we  may 
describe  the  motion  in  complex  notation  as  being  the  real  part  of 


wj(xj»  =  X<xj>  ell0t  (2) 
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Figure  2.  Massless  Beam  Segment 


Substituting  into  Equation  (1),  we  have 


(EDj 


d4  X  (Xj) 

dxj4 


iw  t 


e 


0 


or 


(3) 


Equation  (3)  obviously  has  the  solution 


X  (xj)  =  Axj  +  B 


+  Cx,  +  D 


(4) 


The  physical  meaning  of  the  arbitrary  constants  A,  B,  C,  and  D  is 
easily  seen  to  be 


D  =  X  (0) 

C  =  X'  (0) 

B  =  \  X"  (0)  =  \ 


M(0l 

(EDj 


A 


X'"  (0) 


J021. 

(EDj 


(5) 


Except  for  the  oscillatory  factor  eiwt,  D  is  the  deflection  at  xj  =  0, 

C  is  the  slope  at  xj  =  0,  B  is  1/2  (El)  j  times  the  moment  at  xj  =  0, 
and  A  is  1/6  (El) j  times  the  shear  at  xj  =  0.  We  have  now  essentially 
accomplished  a  first  step  in  the  transfer  matrix  procedure.  We  have 
related  the  deflections  and  internal  force  quantities  which  describe  the 
condition  or  state  of  the  structural  segment  at  a  boundary  point  of  that 
segment  to  the  state  at  any  point  in  that  segment.  In  the  strict  sense, 
the  mechanical  state  of  the  straight  beam  segment  is  completely 
specified  by  the  continuous  deflection  function.  However,  in  this  and 
later  examples,  we  will  when  we  refer  to  the  state  of  an  elastic  element 
mean  not  only  the  necessary,  basic  deflection  components,  but  also  the 
internal  force  and  other  deflection  components,  (or  equivalently,  the 
derivatives  of  the  deflection  components)  necessary  to  calculate  the  basic 
deflections  at  one  point  when  they  are  known  at  another.  We  complete 
this  first  step  by  substituting  Equations  (5)  into  Equation  (4)  and 
differentiating  thrice  to  obtain 
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+  0(0)  X.  +  X(0) 


V(0) 

v5 

M(0) 

x  (Xj) 

"  6  (EI)j 

xj 

2(EI)j 

e(xj) 

.  V(0) 

2(EI)j 

+ 

MiO).  . 
(El). 

M  (x^ 

=  V(0)x  + 

M(0) 

Vfxj) 

=  V(0) 

Equations  (6) 

will  be  the  basis 

of  what  we 

(6) 


matrix. 


To  show  by  example  that  other  sufficient  analyses  of  the  struc¬ 
ture  will  lead  to  these  same  relations,  we  will  now  note  that  Equations  (6) 
could  have  been  derived  from  the  elementary  theory  of  strength  of 
materials.  By  so  doing,  by  means  of  the  following  diagrams,  we  will 
explicitly  state  our  sign  convention.  First,  to  derive  the  shear  relation 


V(0) 


V(xj) 


V(0) 


To  derive  the  moment  relation 


M  (0) 


V(0) 


|M(Xj)  =  M(0)  + 


Xj  V(0) 


Vlxj) 


If  the  moment-area  method  is  used  to  derive  the  expressions  for 
0  (xj)  and  X  (xj)*  we  have 
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Once  again 


b<xj> 


0(0)  + 


M  (0) 

(EI)j 


.  V  (°? 

2(EI)i 


x, 

J 


z 


X<xj) 


X(0)  + 


0  (0)  x. 


+ 


-MjO)  » 

2  (El).  *j 


-_YIP) 

6  (  E  I) . 


If  we  now  specialize  Eqs.  (6)  by  letting  xj  =  1 we  may  view  the 
result  as  a  means  of  carrying  over  or  transferring  the  mechanical 
state  of  the  beam  at  the  right  hand  side  of  station  j-1  to  the  left  hand 
side  of  station  j.  We  may  easily  arrange  this  series  of  transfer 
equations  in  matrix  form  as  follows 


x-w 


0 

M 

V 


= 


J 


3  '~ 


0 

0 

0 


1 

0 

0 


i2 

i3 

2EI 

6EI 

l 

JL 

El 

2EI 

1 

l 

0 

1 

'l 


r 


r  (6*> 


j  j-i 
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where  both  the  superscripts  (i  =  left  of,  r  =  right  of)  and  subscripts  on 
the  column  matrices  (called  state  vectors)  refer  to  a  station,  and  the 
subscript  on  the  square  matrix  indicates  that  the  quantities  that  com¬ 
pose  the  elements  of  this  matrix  are  those  of  the  indicated  beam  seg¬ 
ment.  As  mentioned  previously,  the  amplitudes  X,  0  ,  M  and  V 
are  not  the  only  choice  as  the  four  elements  of  our  state  vector.  (For 
example,  we  could  have  selected  X  and  its  first  three  derivatives.) 
This  choice,  however,  is  much  more  convenient  since  all  of  these 
elements  are  of  interest  in  themselves,  and,  more  importantly,  they 
will  facilitate  the  introduction  of  the  boundary  conditions  into  the 
final  equation.  For  the  sake  of  being  more  concise,  we  rewrite 
Eqs.  (6)  as 

tz>  *  =  [Flj  (Z) '  j  (7) 


where  the  respective  definitions  of  terms  are  obvious.  The  matrix  [  F  ] 
is  called  a  field  transfer  matrix.  We  remark  that  if  the  displacement 
and  force  elements  of  the  state  vector  are  separated  and  the  corres¬ 
ponding  deflection  and  force  elements  are  arranged  in  mirror  image 
position!  and  if  a  suitable  sign  convention  is  chosen,  then  the 
associated  transfer  matrices  are  usually  symmetrical  about  the 
cross-diagonal.  This  is  the  cas  e  here  and  in  the  following  two 
examples.  Another  interesting  feature  of  field  transfer  matrices  is 
that  their  inverses  can  always  be  calculated  by  simply  writing  the 
field  transfer  matrix  for  the  opposite  direction. 


To  have  complete  freedom  of  description!  another  basic  trans¬ 
fer  matrix,  called  a  point  transfer  matrix,  is  needed.  It  transfers  the 
state  vector  across  a  station.  Again  let  us  consider  the  most  general 
case;  i*  e.  ,  let  there  be  a  mass,  a  translational  and  a  torsional  spring 
at  station  j.  The  free-body  diagram  of  station  j  is  as  shown  in 
Figure  3. 


Utilizing  d'Alembert’s  principle,  we  may  write: 


M. 


M.  +  (-0J2  J.)  .0  .  +  k.  e . 
J  J  J  J  J 


v,  -  (-U)2u  XJ  -  k,X 


J  i 
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Figure  3.  Forces  and  Moments  Acting  as 
a  Discrete  Mass 


We  also  may  write,  because  of  the  continuity  of  the  deformation 
across  station  j : 


0 


r 

i 


Again  arranging  our  result  in  matrix  form,  we  have 


r  \ 
X 

r 

0 

0 

0 

e 

*  _ 

0 

1 

0 

0 

M 

0 

k-u2J 

1 

0 

V 

-  > 

j 

-k.j+(jd2U 

0 

0 

1 

M 


(8) 
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or 


(Z)  r  =  [G].  {Z}f  (8a) 

j  j  j 

We  are  now  ready  to  undertake  the  analysis  of  any  undamped 
Myklestad  beam  possessing  any  end  conditions.  To  make  clear  the  gen- 
eral  case,  we  choose  the  specific  beam  shown  in  Figure  4  where  any 
combination  of  elastic  restraints  may  be  present  at  the  three  inter¬ 
mediate  stations. 


Figure  4.  Example  Myklestad  Beam 

Proceeding  from  left  to  right,  we  may  by  successively  applying  Eqs.  (7) 
and  (8)  write 

t Z }l  =  [FI  {Z}  r 
1  1  0 

tz)'  =  [gi1  m-!  (z}tq 

«>1  *  IcIj  •”  to] !  Tb-J  1  (Z)  l 

or  more  concisely 

-  ‘4l  T]0r  {Z )l  (9) 

We  call  [  T]  a  transfer  matrix;  i.e.,  no  adjective.  Its  superscripts 
and  subscripts  indicate  the  extent  of  the  transfer  of  the  state  vector. 

The  boundary  conditions  of  our  example  are: 

at  (i,  4)  X  =  G  =0 

at  (r,  0)  X  =  M  =0 
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(A  similar  result  would  ensue  for  other  boundary  conditions.)  Insert¬ 
ing  these  boundary  conditions  in  the  matrix  equation  (9)  we  have 


< 


M 

V 


4 


C11 

C12 

C13 

fc14 

c2l 

t22 

t23 

c24 

C31 

t32 

t33 

c34 

C41 

c42 

c43 

c44 

r 


0 


< 


0 

0 

0 

V 


r 

>0 


from  which  we  note  that  we  may  extract  two  linear  homogeneous  equa¬ 
tions,  i.  e.  ,  the  submatrix  equation 


The  existence  of  a  nontrivial  solution  to  Eq.  (10)  requires  that  the 
determinant  of  the  coefficients  be  equal  to  zero.  Since  these  coef¬ 
ficients  are  all  functions  of  the  free  vibration  frequency,  our  equation 
for  the  natural  frequencies  of  the  system  is  then 


i 


M  2 


"■l  4 


r 


4 


*2  2 


* 2  4 


0 


=  A  (go  )  =  0 


(11) 


Since  A(w)  will  for  this  model,  be  a  polynomial  in  to 2  of  order  equal 
to  the  number  of  intertia  parameters,  both  rotary  and  translational,  a 
numerical  rather  than  an  exact  solution  for  w  is  to  be  obtained  in 
almost  all  cases.  Of  course,  all  the  roots  for  go2  are  real  and  positive. 
Once  the  natural  frequencies  are  known  from  equation  (11),  it  is  a  simple 
matter  to  compute  the  force  and  deflection  mode  shapes  for  this  beam. 

We  start  with  Eq.  (10)  which  tells  us 


(12) 
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If,  for  example,  we  normalize  (i.  e.  assign  a  unit  value  to)  V,  we  know 
all  four  components  of  the  modal  state  vector  at  (r,  0).  It  is  now  only 
necessary  to  use  the  transfer  matrices  [  G (oo  n)  ]  and  [F(xj)]  to  com¬ 
pute  the  modal  state  vector  at  any  point  of  continuity  on  the  beam. 

The  above  steps  are  the  basic  elements  of  the  transfer  matrix 
approach.  In  review,  we  start  at  one  point  of  the  structure  and 
proceed  to  an  adjacent  point  using  whatever  state  vector  is  necessary 
to  describe  plus  transfer  the  mechanical  state  of  the  structure.  The 
transfer  matrices  that  carry  the  state  vector  from  one  point  to  another 
can  be  calculated  by  solving  the  overall  equation(s)  of  motion  of  each 
structural  segment.  When  the  proper  field  and  point  transfer  matrices 
are  determined,  a  transfer  matrix  from  one  boundary  to  another  can 
be  constructed,  and  then  partitioned  to  solve  for  the  natural  frequencies  and 
mode  shapes. 

The  present  analysis  could  be  expanded,  for  instance,  to  con¬ 
sideration  of  the  vertical  vibrations  of  a  cantilevered  beam  with  a 
horizontal  bend.  We  could  further  slightly  complicate  the  matter  by 
placing  the  concentrated  masses  and  rotary  inertias  on  rigid  extensions 
from  the  beam  axis.  This  model  might  be  a  suitable  model  for  a 
cantilevered  high  aspect  ratio  wing.  See  Figure  5. 


Figure  5.  Bent  Beam 
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We  won't  stop  to  analyze  this  model,  but  only  comment  on  those 
of  its  features  that  are  not  found  in  the  previous  example.  First, 
simply  because  of  the  horizontal  bend,  a  vertical  force  or  bending 
moment  on  the  outboard  section  of  the  beam  will  produce  a  torsional 
moment  (and  deflection)  as  well  as  a  bending  moment  on  the  inboard 
section.  Therefore  we  would  have  to  include  the  local  beam  twist  and 
torsional  moment  in  an  analysis  of  this  structure.  Secondly,  even 
without  the  horizontal  bend,  we  would  be  compelled  to  consider  tor¬ 
sional  twist  and  torsional  moment  because  the  vertical  accelerations 
of  the  offset  concentrated  masses  would  produce  inertial  torques 
about  the  beam  axis.  Therefore,  our  model  would  properly  include 
not  only  the  rotary  inertias  associated  with  bending  slope  deflec¬ 
tions,  but  also  the  rotary  inertias  corresponding  to  torsional  de¬ 
flections,  and  cross  products  of  inertia. 

We  could  start  our  analysis  at  the  beam  tip  and  work  our  way 
toward  the  bend  just  as  in  the  previous  example.  (The  additional  equa¬ 
tion  of  motion  for  a  massless  beam  segment  is:  the  effective  torsional 
stiffness  times  the  second  derivative  of  the  twisting  angle  equals  zero. 
The  equations  of  motion  of  beam  stations  can  again  be  derived  using 
D' Alembert’s  principle.)  At  the  bend  we  simply  connect  the  locally 
oriented  state  vectors  on  each  side  of  the  bend  by  a  matrix  that  resolves 
the  outboard  set  of  deflections  and  internal  forces  into  components 
referenced  to  the  inboard  coordinate  system.  (Matrices  are  a  rather 
efficient  way  of  handling  coordinate  rotations.) 

We  are  not  limited  to  models  where  the  inertia  properties  are 
made  discrete.  Let  us  now  investigate  how  we  may  apply  the  transfer 
matrix  method  to  a  beam  with  distributed  inertia  properties. 
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in.  THE  CONTINUOUS  BEAM  PROBLEM 


Consider  the  case  of  a  uniform  beam  segment  of  symmetric  cross 
section  and  of  uniform  density  between  two  stations,  say  j~l  snd  j . 
See  Figure  6. 


f 


Figure  6.  Continuous  Mass  Beam 


The  equation  of  motion  of  this  segment  can  be  derived  from  the  basic 
beam  bending  equation.  It  is,  for  vibrating  motion 

El  -T—r  =  -  pAw  (13) 

9  x 

where  p  is  the  mass  density,  A  the  cross-sectional  area,  and  each 
dot  above  a  symbol  indicates  one  partial  differentiation  with  respect 
to  time.  We  have  dropped  the  subscript  j  from  our  various  quan¬ 
tities  for  the  sake  of  simplifying  the  writing  of  our  equations.  If  we 
wish  to  examine  a  single  frequency  sinusoidally  varying  motion,  we 
may  again  describe  the  vertical  deflection  as 


w{x,  t)  =  X  (x)  e 


Thus  our  equation  of  motion  reduces  to 


El 


d4  X 
dx4 


=  p  A  w2  X 


This  equation  may  be  rewritten  as 


d4X 

d^ 


ir4  X  =  0 


(14) 
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(14a) 


where  o4  =  P.Aa)2 

El 

Four  independent  solutions  to  the  differential  equation  (14)  are  exp(ox), 
exp(-ox),  exp(iox),  and  exp(-iax).  Note  that  these  solutions  are 
independent  regardless  of  which  of  the  four  roots  of  a4  is  used,  and  in 
fact,  the  same  four  functions  are  obtained  regardless  of  which  of  the 
four  roots  is  used.  For  later  convenience  we  will  specify  the  use  of 
the  positive  real  root.  Now  any  four  independent  linear  combinations  of 
the  above  four  solutions  are  also  a  complete  set  of  solutions.  The  most 
advantageous  set  of  combinations  is 


(x)  = 

1 

2 

(cosh  <rx 

+ 

cos  <rx) 

(x)  = 

1 

2 

(cosh  <rx 

- 

cos  <rx) 

(x)  = 

1 

2 

(sinh  ax 

+ 

sin  or x ) 

(X)  = 

1 

2 

(sinh  <rx 

sin  crx) 

The  advantages  of  this  set  of  combinations,  the  reasons  for  their 
selection,  are 

(1)  f0  and  f2  are  real  even  functions  of  x 
fj  and  f3  are  real  odd  functions  of  x 

(2)  The  index  order  rotates  under  differentiation  or  integration;  e.  g. 

V  (x)  =  cr  f3  (x) 

f,  '  (x)  =  cr  f0  (x) 

f2  '  (x)  =  <r  f,  (x) 

f3'  (x)  =  <rf2  (x) 

and  (3)  f  0  (0)  =  1;  f^  (0)  =  0  for  j  4  0. 

Then,  when  we  write  our  solution  in  the  standard  form 

X(x)  =  A0  4  (x)  +  Ai  *i  (x)  +  A2  f2  (x)  +  A 3  f3  (x)  (16) 

we  can  see  by  virtue  of  items  (2)  and  (3)  above  that  we  obtain  the  happy 
result 


16 


Ao  =  X(0) 
o-  A,  =  X'(0) 

or2A2  =  X"(0) 
cr3  A3  =  X'n{0) 


or 


\  = 


0(0) 


or 


or 


A  _  JIM 


i r2  El 


A.  = 


y(Q) 

rEI 


Therefore 

0  (0) 

X(x)  *  X(0)  f0  (x)  +  — —  f,  (x) 


+ 


M 

a*  El 


U  (x)  + 


ViO} 

rEI 


U  (X) 


(17) 


(l6a) 


So,  once  again  by  obtaining  a  suitable  solution  to  the  equation  of  motion 
of  an  elastic  element,  we  have  obtained  the  basic  relation  for  the  field 
transfer  matrix  of  that  element.  The  remaining  relations  are,  of 
course,  obtained  by  differentiating  Eq.  (l6a) 


X'  (x)  =  <rX(0)  f3  (x)  + 

+  2US>  fl  (X)  + 

o-  EX  1  y  ' 


©  (0)  f0  (x) 

1‘  <x> 


X"  (x)  =  cr2  X  (0)  f2  (x)  +  0-0  (0)  f3  (x) 

+  ^e^  f°  (x)  +  Vm  f>  (x) 


X'"(x)  =  (T3  X  (0)  f,  (x)  +  .r2  ©  (0)  f2  (x) 

via 


M  (0)  ,  ,  % 

+  O'  -eP  f3  W  +  EI 


fo(x) 
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Combining  these  results  for  x-  =  I  ■  in  matrix  form,  we  have 
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(18) 


j-1 


Note  that  this  field  transfer  matrix  is  also  cross-symmetrical. 


At  this  point,  let  us  include  damping  effects  as  part  of  our 
analysis.  We  do  so  in  the  realization  that  for  a  damped  system 
sinusoidal  motion  exists  only  when  the  system  is  subjected  to  a 
sinusoidal  excitation.  Structural  damping  at  joints  and  in  elastic 
segments  can  easily  be  included  in  the  conventional  manner  in  our 
equations  of  motion.  At  joints,  that  is  at  the  location  of  elastic 
springs,  we  introduce  structural  damping  by  simply  replacing 
k  by  k(l  +  16)  and  k  by  k  ( 1  +  i  c)  where  6  and  £  are  structural 
damping  factors.  Similarly,  structural  damping  in  an  elastic 
element  is  described  by  replacing  El  by  EI(1  +  iy),  where  y  is 
another  structural  damping  factor.  Of  course,  these  changes  in  the 
equations  of  motion  are  carried  through  in  identical  form  to  the  point 
and  field  transfer  matrices.  We  may  also  take  into  account  an 
equivalent  viscous  damping  by  adding  another  term  to  our  equations 
of  motion.  For  example,  the  equation  for  an  elastic  segment  becomes 

EI(1  +  iy)  +  cw  +  p  Aw  =  0 

with  no  external  excitation  on  this  particular  segmeit.  To  solve  this 
equation  we  again  let  w  =  X  (x)  exp  (ioo  t) ,  and  substitute  to  obtain 


d4  X  p  w2  A  —  lew 

dx4  “  E I  ( 1  +  i  y) 


(19) 
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We  can  return  to  the  form  of  equation  (14)  by  simply  redefining  &  as 


04  =  pto2  A  -  icfa)  (19a) 

El (1+iy) 

Thus,  again,  the  field  transfer  matrix,  Eq.  (18)  retains  the  same  form. 
In  this  case  all  the  quart ic  roots  are  complex.  As  before  it  is  permis¬ 
sible  to  use  any  one  of  the  four  roots  as  o  in  the  computation  since  we 
only  require  that  sinhox,  coshox,  sinox,  and  cosax  be  independent  of  each 
other. 


Thus,  since  the  form  of  the  transfer  matrices  remains  unaltered 
by  the  presence  of  damping  effects,  in  our  subsequent  examples  we  will 
take  for  granted  the  presence  of  the  various  damping  parameters  in 
their  appropriate  locations. 

As  a  final  comment,  before  we  pass  to  the  topics  of  forced 
vibration  and  skin-stringer  construction,  the  method  of  transfer  matrices 
can  also  be  applied  to  beam  grids  and  frames [2].  Very  briefly,  a  single 
path  through  the  structure  is  chosen  and  by  means  of  the  continuity  and 
equilibrium  equations  of  the  joints  and  the  boundary  conditions  at  beam 
ends  other  than  the  two  that  mark  the  beginning  and  end  of  the  chosen 
path,  the  structure  that  lies  off  the  chosen  path  is  represented  by 
point  transfer  matrices  on  the  chosen  path. 
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iv.  response  to  forced  vibration 


Before  we  see  how  to  use  transfer  matrices  to  calculate  the 
response  of  a  structure  to  forced  vibration,  let  us  pause  to  review 
a  bit  of  the  basic  theory  of  mechanical  vibrations. 

The  response  to  a  structure  under  forced  vibration,  either 
deterministic  or  random,  can  be  described  by  either  of  the  following 
two  basic  functions  [6]; 


the  impulse  response  function 
the  frequency  response  function 


h 

pq 


(t) 

(w) 


The  respective  meanings  of  these  two  functions  are  as  follows.  Let  the 
structure  be  at  rest  at  t  =  0,  and  let  the  excitation  be  a  unit  impulse* 
applied  at  station  q  at  time  t  =  0.  Then  hpq  (t)  is  the  response  at 
station  p.  See  Figure  7  below.  Let  the  structure  be  stable  (say 
positively  damped),  and  let  the  excitation  be  a  unit  sinusoidally  varying 
load  (e1W  )  applied  at  station  q.  Then  the  steady  state  response  at 
station  p  is  Hpq  (w)  eiwt. 


Figure  7*  Beam  of  General  Configuration 


*  A  unit  impulse  at  time  x=t  is  described  mathematically  by  the 
Dirac  delta  function,  Its  definition  is,  for  an  arbitrary 

function  F(t) 

ft+C 

r c t )  =  r c t ) *  6(t  -  T>dx  =  r(T>*  6(t  -  x)dx 

J_co 

where  e  is  arbitrarily  small.  It  is  sometimes  described  as  a  func¬ 
tional  because  it  is  not  a  function  in  the  ordinary  sense,  i.e.  while 
<$(t)  -  0  for  almost  all  t,  it  is  not  specifically  defined  for  t  within 
the  interval  (t-e,  t  +s) ,  e  arbitrarily  small. 
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An  arbitrary  forcing  function  can  be  considered  as  a  continuous 
sequence  of  impulses,  i.  e.  we  may  write 

P  (t)  =  J+°°  P  (t)  6(t  -  T)dT 
q  •'-cd  q 

A  response  at  station  p  to  the  forcing  function  can  be  described  with 
the  aid  of  the  appropriate  impulse  response  function  for  the  two  stations, 
i.  e. 

°p  (t)  =  P,  (T)  hpq  (t  -  T)dT 

On  the  other  hand,  if  a  spectral  analysis  of  the  forcing  function  is  per- 
formed  so  that 

_  ,  .  -+oo  —  t  .  iwt  , 

p  (t)  =  /  P  (w)  e  du 

q  -oo  q 

where 

P  («)  =  -j-  /+°°  P  (t)  e"1Wt  dt 
q  '  2tt  _oo  q 

then  the  response  to  the  forcing  function  Pq  (t)  can  be  written 

fi  (t)  =  f°°  Pn  (w)  H  (co)  elwtdw 
P  '  J-oo  <1  pq 


Using  a  few  simple  steps  we  can  show  that 


H  (w)  =  /°°  h{t)  e"lliJt  dt 

pq  -oo 


h  (t)  =  J°°  H  (w)  e^dw 

pq  2tt  j_oo  pq 


That  is, outside  of  the  adjustment  of  a  constant  factor,  these  functions  are 
Fourier  transform  pairs. 

A  simple  example  of  these  functions  is  to  be  found  in  the  case  of 
a  single  degree  of  freedom  system.  Here  p  =  q  =  1  (omitted),  and 
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h(t)  = 


-  i 


exp 


mu  /l-C  *  '  '  n  n 
n 


(-C«  t  +  iu  A  -  42  t),  for  t  >  0 


0,  for  t  <  0 


and 
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m  (w  2  -  oj2  +  2  i  £,  w  to  ) 

n  *  n' 


The  technique  of  transfer  matrices  can  be  used  to  compute 
frequency  response  functions.  For  illustration,  consider  the  non- 
uniform  beam  with  end  conditions  sketched  below. 


Figure  8.  Segmented  Beam  of  General  Configuration. 


The  beam  is  approximated  by  a  finite  sequence  of  uniform  sections 
between  stations.  Any  arrangement  of  interior  elastic  supports  may 
be  included  in  the  problem.  We  are  interested  in  obtaining  the  response 
at  station  k  due  to  a  unit  sinusoidal  force  at  station  j.  We  may  need 
to  be  more  specific  about  the  response,  that  is,  decide  whether  the 
response  on  the  left  hand  side  or  the  right  hand  side  of  station  k  is  the 
response  to  be  calculated.  For  the  sake  of  definiteness,  let  us  decide 
we  want  to  know  the  response  on  the  right  hand  side. 
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To  develop  the  equations  in  a  natural  way,  let  us  again  start  at 
the  left  end  of  the  beam  and  proceed  to  the  right.  By  our  previous 
work,  we  can  immediately  work  our  way  to  the  left  hand  side  of  the 
impressed  force. 


(Z> 


3 


t  r 
[  T  ] 


j 


{Z} 


The  impressed  force  is  simply  introduced  by  writing 


(Z) 


j 


r  ,  r  r 
[  T ]  {Z}  + 


0 

KiJ 


{Recall  the  factor  eiwt  is  common  to  all  terms  of  these  equations.) 
Proceeding  from  this  point,  we  have 


so  that 


r  r  r  r 

(Z)  =  [  T  ]  {Z> 

k  k  j  j 


r  r  r  r  r,  .  r 

{ Z )  =  [  T  ]  (Z)  +  [  T  ] 

k  k  o  o  k  j 


/*0~\ 
0 
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(20) 


This  is  our  first  of  two  basic  equations.  In  expanded  form  it  is 
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(20a) 


where  we  have  inserted  the  boundary  conditions  at  station  0.  In 
addition  to  the  elements  of  the  state  vector  at  (r,  k),  the  slope  and 
shear  at  (r,  o)  are  unknown.  These  elements  can  be  eliminated  by  use 
of  the  boundary  conditions  at  station  N  in  our  second  basic  equation, 
which  is  merely  a  continuation  of  the  first  equation.  It  is 
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(21) 


Now,  from  the  above  two  (4  x  1)  matrix  equations,  (20a)  and  (21),  we  can 
extract  the  following  two  equations,  respectively. 


Then  solving  the  second  of  these  for 


and  substituting  this  quantity  into  the  first,  we  obtain 
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< 


M 

V 
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*2  2 
*3  2 
*4  2 


k 


tl4 
*2  4 
*3  4 
*4  4 


(22) 


The  left  hand  side  is  the  set  of  amplitudes  of  the  response  as  functions 
of  the  frequency  w  due  to  a  unit  sinusoidal  force  excitation,  which  by 
definition  are  the  frequency  response  functions  for  these  circum¬ 
stances.  Of  course,  the  frequency  response  functions  for  a  unit 
sinusoidal  moment  can  be  obtained  in  an  exactly  analogous  way.  The 
result  is 


*12 
*2  2 
*3  2 
*4  2 


(22a) 


Of  course,  other  boundary  conditions  lead  to  the  selection  of  other 
elements  of  the  transfer  matrices. 
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v.  FLAT  STRINGER -PANEL  SYSTEMS 


We  will  discuss  two  types  of  stringer-panel  systems.  We  will 
begin  with  flat  stringer -panel  systems,  and  when  we  have  determined  the 
method  of  solution  for  these,  we  will  proceed  to  discuss  the  more 
eonplex  case  of  curved  panel  systems.  A  sketch  of  a  typical  flat  panel 
array  is  shown  in  Figure  9* 


7 

b 

I 


rm 

0  12 
^  Stringer  No. 


n  '  nT’HL 


Figure  9.  Flat  stringer-panel  model 


To  simplify  the  presentation  of  both  the  flat  and  curved  panel 
analyses,  we  will  limit  our  attention  to  free,  undamped  vibrations. 
The  treatment  of  the  cases  of  damped  and  forced  vibrations,  as  has 
been  explained,  is  little  different.  A  simple  forced  vibration  example 
will  be  given  in  Section  XII. 
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Note  that,  while  in  the  case  of  beams  we  made  no  restric¬ 
tions  concerning  the  boundary  conditions,  in  our  flat  panel  mathe¬ 
matical  model  we  specify  simple  support  boundary  conditions  at  the 
frames.  As  will  be  seen,  this  particular  arrangement  allows  us  to 
reduce  this  essentially  two  dimensional  problem  to  a  one  dimensional 
problem.  As  pointed  out  before,  once  we  have  a  one  dimensional 
problem,  we  may  employ  the  method  of  transfer  matrices.  Most 
actual  construction  is  not  such  that  the  conditions  of  simple  sup¬ 
port  exist  at  the  frames.  However,  because  the  distance  between 
stringers  is  usually  one  half  or  less  than  one  half  the  distance 
between  frames,  and  since  fatigue  failure  can  be  expected  to  begin 
in  the  panel  skin  adjacent  to  the  center  of  a  stringer,  this  model 
is  quite  practical  for  the  study  of  the  important  problem  of  fatigue 
stresses  in  actual  panel  construction.  More  generally,  this  same 
argument  can  be  extended  to  claim  validity  for  applying  the  results 
of  this  model  to  construction  with  different  boundary  conditions  at 
the  frames  if  the  response  of  interest  is  located  near  the  mid- 
distance  line  between  frames. 

The  standard  plate  bending  equation  of  motion  for  the  panel 

skin  is 


34w 

3x4  +  Bx2  By2 


+ 


(23) 


For  a  discussion  of  the  underlying  assumptions  and  derivation  of 
Eq,  (23)  see  Reference  [7],  pages  39-1  through  39-3,  In  this  problem, 
with  simple  support  at  the  frames,  we  may  obtain  a  solution  in  the 
form 


w(x,  y,  t)  =  Yn  (y)  sin  BI* 


(24) 


where  we  take  advantage  of  the  known  modal  form  in  the  x-direction. 

For  each  n,  the  function  Yn  is  the  amplitude  function  in  the  y-dlrection 
It  corresponds  to  X(x)  in  the  beam  problem,  and  this  unknown  is,  of 
course,  only  a  function  of  one  space  coordinate. 

Substituting  Eq,  (24)  into  Eq,  (23),  after  cancelling  exp(iwt) 
and  sin(mTrx/b),  we  obtain  for  each  stringerwise  modal  number  n 
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(25) 


Y 


n 


We  have  two  unknowns  in  the  one  equation,  w  and  Yn.  We  need 
to  determine  those  special  values  of  co  (in  mathematical  terms  the 
eigenvalues,  or  in  dynamical  terms  the  natural  frequencies)  for 
which  we  can  also  determine  solutions  for  Yn  up  to  an  arbitrary  multi¬ 
plicative  constant,  (the  eigenfunctions,  or  normal  modes  of  an  undamped 
system. ) 


We  can  write  the  characteristic  roots  of  Eq.  (25)  as  +<Tj  ,  -<Tj  , 


(26) 


2 

Since  for  each  n,  (Dl/4mr/hl/4  p  1/4  b)  is  the  smallest  natural 
frequency  of  an  infinitely  long  panel  row  without  stringers  L  8 -I  ,  and 
because  the  addition  of  stringers  and  boundary  conditions  at  a  finite 


distance  would  increase  the  stiffness  and  hence  the  natural  frequen¬ 
cies  of  such  a  built-up,  finite  system  ,  we  can  conclude  that  <*2  is  a 
real  quantity. 


Let  us  now  turn  to  our  transfer  matrix  techniques.  We  are 
here  interested  in  determining  the  deflection,  slope,  moment,  and 
shear  in  the  panel  skin  along  any  plane  perpendicular  to  the  x-axis. 
The  deflection  and  slope,  of  course,  depend  upon  Yn  and  Y'n 
respectively.  Sincet^J 


D( 


3  *  w 


V  =  +  D 

y 


3 

3  y 


V 


+ 


32w  N 

T?) 


32w\ 

T?  ) 


27 


we  see  that  we  can  describe  these  quantities  of  interest  by  employing 
Yn  and  its  first  three  derivatives.  Therefore,  as  a  stepping-stone  to 

fYUr  v»Mte  V®c^?r  LYn  Yn  Mn  Vn-I>  let  us  consider  the  state  vector 

L:n  Yn  Yn  Ynll  anc*  discover  how  it  may  be  transferred  across  a  field 
i.e.  across  a  panel. 


We  already  have  the  first  part  of  the  solution  to  the  panel  prob- 
lem  in  terms  of  the  solutions  for  the  chatacteristic  roots  of  our  reduced 
equation  of  motion,  Eq.  (26).  These  roots,  of  course,  allow  us  to 
write  a  solution  for  Yn  in  terms  of  (e.g.  expodential  functions  of)  o> 
for  each  panel.  Then,  analogous  to  the  continuous  beam  problem,  our 
next  step  is  to  express  the  panel  solution  in  a  form  appropriate  to  our 
state  vector  |_Yn  Yn  Yn  Yn_|  .  Again  we  start  by  expressing  our  solution 
in  the  form  of  the  sum  of  four  independent  functions,  i.e. 


'n(Y)  =  Ao  fo<y>  +  h  fi(y)  +  A2  f2(y)  +  a3  f3(y) 


(27) 


so  that  among  other  conveniences,  A„  =  Yn(0),  A,  =  Y'fOl  A-.  =  Y'^ni 
and  A3  =  Yn"’(0).  As  in  the  cas^  of°the  Seam,’  we \y-pass  solStionS  [n  ’ 
the  form  of  exponential  functions,  and  concentrate  on  combinations  of 
the  hyperbolic  and  circular  functions.  Specifically,  these  independent 
soiuttons  are  sinh  Ojy,  cosh  ojy,  sin  o2y,  and  cos  o2y.  The  combination 
of  these  that  we  seek  is 


f0(y) 


°i2  +  °22  2 


(o9^  cosh  Oiy  +  a? 


°iy  +  ai  cos  °2  y) 


fi(y)  = 

f2(y)  = 
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Slnh  °iy  +  57 


O-i  +  Or 


(cosh  o^y  -  cos  02y) 
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02  y) 


(28) 
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2  0  2  (  0  Slnh  °iy  "  —  sino2y) 


(A  general  procedure  by  which  we  can  calculate  the  desired  functional 
forms  of  our  solution  will  be  illustrated  for  the  more  complex  case  of 
a  curved  panel  segment.)  The  functions  f  and  f2  are  even  while  f, 
and  f3  are  odd,  fs(r)(0)  =  0  when  r  i  s,  and  fs(r)  =  l  when  r  =  s,1 
for  r,  s  -  0,1,  2,  3.  However,  because  a-,  ^  02,  ff(y)  ^  f  .  1  (y) ; 

i-e*  the  nice  index  rotating  property  that  was  present  in  the  case  of 
the  beam  problem  does  not  exist  here.  As  a  result  of  this,  we  will 
requirevmore  than  four  functions  to  construct  the  transfer  matrix  for 

n  J* 
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Now  that  we  have  as  our  solution  Eqs.  (28),  we  can  write 


Yn  (y)  =  Yn  (0)  f0  (y)  +  (0)  f,  (y)  +  Y^  (0)  f2  (y) 

+  Y»«  (0)  f3  (y) 


(27a) 


Thus,  again,  by  differentiating  and  then  setting  y  =  X^,  we  can  establish 
the  relations  between  Yn  and  its  derivatives  at  y  =  X  .  on  one  hand, 
and  Yn  (0)  ,  Y^  (0),  Y|{  (0),  and  YJJ'  (0)  on  the  other  hand.  These 
relations  in  matrix  form  are 
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where  y  =  X.  corresponds  to  (f,  j)  and  y  =  0  corresponds  to  (r,  j-1), 
and  J 
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where 


C0  =  fb  (X.) 
S.,=  f,  (X.) 
C_2=  f2(x.) 

S-3  =  (X.) 
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At  this  point  we  will  convert  the  state  vector  |_Yn  Y”  Y^'J 
to  one  which  will  represent  the  variations  in  moment  and  shear  in  the 
y  direction.  We  do  this  by  expressing  the  moment  and  ~u" --  in  fhp 
same  form  as  the  vertical  deflection,  i.e. 


M  =  eiwt  Mr,  (y)  sin  BIB. 

y  b 

V  =  eiWt  Vn  (y)  sin  BIB. 

y  n  .b 


(30) 


From  Eqs.  (27)  and  (30),  by  means  of  the  orthogonality  properties  of  the 
sine  function,  we  can  obtain 

MnW  "  D  [Yn  *  ’'(f)  ‘  V 

(31) 

V  (y)  =  D  [  Y"‘  -  (2  -  v)(^  )2  Y'  ] 
n  n  \  b  J  n 

This  immediately  leads  to  the  matrix  form 
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where  the  conversion  matrix  and  its  inverse  in  detail  are 
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Now  we  are  in  a  position  to  write 
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where  [  F  1  •  =  [b]j[r]j  [  B  ]  j  is,  of  course,  the  field  matrix 
which  transfers  the  state  vector  LYn  Y^  Mn  V^J  across  a  panel. 
The  matrix  [  F  ]  is  symmetrical  about  its  cross-diagonal,  and 
[  F(bn)  ]  _1  =  [  F  (-bn)  ]  . 

Now  that  we  can  leap  panels  as  we  wish,  we  turn  to  our  last 
hurdle,  the  crossing  of  stringers.  To  find  the  point  transfer  matrix 
which  transfers  a  state  vector  across  a  stringer,  we  must  take  a 
closer  look  at  the  stringer.  The  stringer  does  not  interfere  with 
the  continuity  of  the  deflections  and  slopes  in  the  skin  on  either  side 
of  the  line  of  attachment  between  the  stringer  and  skin.  Thus,  very 
simply, 


Yn  (r,  j)  =  Yn  (f,  j) 

Y;  (r,  j)  =  Y^  (f ,  j) 

However,  the  stringer,  because  of  its  elastic  and  inertial  properties, 
does  produce  what  we  will  consider  to  be  an  abrupt  change  in  the 
moment  and  shear  in  the  skin  at  the  line  of  attachment.  To  determine 
the  form  of  the  change  in  the  shear  and  moment  component*  of  the 
state  vector  we  will  make  use  of  an  existing  theory  which  deals  with  the 
bending  and  torsion  of  a  general  thin-walled  member  of  open  cross- 
section.  This  theory  is  the  result  of  a  number  of  distinguished 
engineers,  a  few  of  which  are:  Timoshenko  (in  1908);  Wagner  (in  1929); 
Goodier  (in  1941);  and  Argyris  (in  the  early  1950's). 

The  basic  deflection  equations  of  the  stringer  are 
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(34) 

EC  4; (4) 
w 

-  GC4j" 

=  r(x) 

where  I  ,  I  ,  and  I  are  the  centroidal  moments  of  inertia  and  pro¬ 
duct  moment3 of  inertia;  G  is  the  shear  modulus  of  elasticity;  the  sub¬ 
script  o  refers  the  deflections  w  and  v  to  the  shear  center  ;  C  is  the 
warping  constant  of  the  stringer  cross  section  with  respect  to  the  shear 
center  ;  C  is  the  Saint- Venant  constant  of  uniform  torsion;  p  and 
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are  the  distributed  loads  at  the  shear  center  along  the  stringer  in  the 
z  and  y  directions  respectively;  7  is  the  distributed  torque  about  the 
shear  center;  and  primes  indicate  differentiation  with  respect  to  x. 

The  first  two  of  equations  (34)  can  be  derived  by  simply  (a)  super¬ 
imposing  the  strains  of  a  stringer  cross-section  due  to  curvatures 
in  the  two  perpendicular  planes  (plane  sections  assumed  to  remain 
plane);  (b)  integrating  the  corresponding  stresses  to  obtain  the  values 
of  the  moments  in  terms  of  the  deflections;  and  (c)  differentiating  to 
obtain  the  distributed  loads.  The  derivation  of  the  third  of  Eqs.  (34) 
can  be  found  in  Reference  [7]  ,  section  36.  8. 

The  distributed  forces  and  the  distributed  torque  about  the  shear 
center  are  a  result  of  the  acceleration  of  the  stringer,  and  (reversing 
our  original  point  of  view)  the  differences  between  the  forces  and 
moments  in  the  skin  on  opposite  sides  of  the  point  of  attachment  S. 
Specifically 


p  (x)  =  (-  Vr  +  V*  )  -pA  w 

c 

q  (x)  =  (N**  -  N*)  -PA  v 
9  9  c 

7  (x)  =  (Mr  -  Mi)  -  (N*-  -  N*)  s 

9  9  z 


(35) 


pAVc  cz 


+  p  A  w  c  —  p  J  *Jj* 
c  y  c 


where  is  the  amplitude  of  circumferential  tension  per  unit  length; 

p  is  the  stringer  mass  density;  A  is  the  stringer  cross  sectional 
area;  Jc  is  the  stringer  area  polar  moment  of  inertia  about  the  cen¬ 
troid  C;  c^,  Cy,  and  &z  are  distances  illustrated  in  Figure  10,  and 
the  subscript  c  on  the  deflections  refers  them  to  the  centroid  C. 

Our  basic  equations  are  expressed  in  terms  of  the  deflections 
at  the  stringer  shear  center  0  and  the  stringer  centroid  C.  Our 
purpose  is,  however,  to  express  ourselves  in  terms  of  the  y-coordinate 
continuous  deflections  of  the  panel  skin.  This  can  be  accomplished  by 
noting  that  the  deflections  of  the  stringer  and  those  of  the  skin  are  the 
same  at  the  point  of  attachment  S.  As  can  be  seen  from  Figure  10, 
the  deflections  of  the  stringer  at  point  S  (no  subscript)  are  geometrically 
related  to  those  at  the  centroid  and  shear  center  as  follows 
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w  =  w  -  s  (1  -  cos  4j  )  »  w 

O  Z  ^ 

v  =  vo  -  8Z  + 

Wc  =  WD  +  Cy  + 

V  =  V  -  C  4* 

c  o  z 


(36) 


Substituting  these  relations  into  our  equations  of  motion  (34)  and  (35),  and 
combining  and  rearranging  them,  we  obtain 


M  - 


i  —  —  .  (4)  r*  n  .i. "  j.  in  s  +  El  s  w(4) 


M  *  EC  di'  -  GC  i(i  +  El  s  v 
ws  S 


Tl4  z 


+  pj  Oj  -  PA  (c  -  s  )  v  +  p  A  cv  w 

a  T  Z  Z  y 


Vi  =  -  El  w(4)  -  El  v(4)  -  El  s  vj<(4) 

T|  Tit,  Tit,  Z 


(37) 


-  p  A  w  —  pAc  ^ 

y 


N*  -  N*  =  El  v(4)  +  El  w(4)  +  El  s  i|>(4) 
<|)  <j>  t,  t,  z 


+  p  A  V  +  p  A  (c  -  s  )  ^ 


z  z 


where  C  =  C  +  I  s  2  is  the  warping  constant  of  the  stringer 
cross- section  w$h  respect  to  S  as  the  center  of  twist,  and  Jg  =  Jc  + 

A  c  2  +  A  (c  -  sz)2  is  the  area  polar  moment  of  inertia  of  the 

stringer  cross-section  with  respect  to  S. 

We  now  recognize  that  the  flat  panel  system  undergoes  negligible 
lateral  motion,  i.  e.  mathematically  v  =  0.  This  is  in  keeping  with  the 
underlying  assumptions  of  our  equation  of  skin  motion,  Eq.  (23).  That 
is,  Eq.  (23)  is  a  bending  theory  equation  in  which  bending  of  the  plate 
produces  no  horizontal  motion  of  points  in  the  plate  middle  plane.  (On 
the  other  hand  a  bending  theory  of  curved  panels  will  obviously  have  to 
include  lateral  deflections.)  So,  substituting  v  =  0  and  ^  =  4*  into 

Eqs.  (37)  we  obtain 
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Mr  -  M*  =  EC _ JLiL_  -  GC 


3w 


WS  3x^y 


3x23y 


+  EInC  s  ^.  +  PI 
TIC  Z  3x4  S 


3 

3y3t5 


+  PA  c 


8fw 
>  3t2 


Vr  -  V8-  =  -  EIn  -  EIn,  b 

n  3X^  n?  z  3x43y 


-  PA  Ifj.  -  pA  cy 


3t' 


5y9 


(38) 


We  conclude  the  derivation  of  the  point  transfer  matrix  by  substituting 
the  expressions  for  w,  M,  and  V,  Eqs.  (24)  and  (31),  into  Eqs.  (38)  to 
arrive  at 


(Mn>r  -  (V*  =  [ECWS  (^)-  +  GC(^)2 
+  [EIne  sz{T)k  -“2PAcyJYn 

=  d  Yn  +  3  Y-n  (39) 

<Vr  -  ~  u2pA]Yn 

-  [EInc  s2  -  w2pA  cy]Yn’ 
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-  dY' 


n 


We  note  that  the  quantities  c,  d,  and  e  are  functions  of  the  natural 
frequency.  Thus,  in  conclusion,  we  can  write 


(z}3  *  l0lJ  <Z,J 

where 


(1*0) 


[G] 


A 


1 


0  0 


0 


0  10  0 


d  c  1  0 


0  1 


A 


and  the  subscript  j  associated  with  [g}  indicates  that  the  elements  of  the 
point  transfer  matrix  are  computed  from  the  physical  data  of  the  j**1 
stringer. 


The  only  exceptions  to  the  cross -symmetry  of  the  stringer  transfer 
matrix  are  the  elements  d  and  -d.  These  exceptions  are  a  result  of  the 
stringer  cross<-section  being  unsymmetrical;  ie.  d=0  if  the  section  is 
symmetrical,  and  when  the  direction  of  transfer  of  the  state  vector  is 
reversed  causing  a  far-side  eccentricity  to  become  a  near-side  eccen¬ 
tricity,  d  becomes  -d. 


We  now  can  proceed  to  calculate  the  natural  frequencies  and  modes 
of  the  entire  panel  array.  Let  [t],  with  proper  subscripts  and  super¬ 
scripts,  denote  the  general  transfer  matrix  relating  the  state  vectors  at 
any  two  stations  on  a  panel  system.  For  example 


and 


(F]k  [0]w  in*.! 


k[T]$  =  [G]k  [F]k  [G]k_i  •••  [F]J+1  [G]j 
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For  the  purposes  of  discussion,  let  us  assume  that  the  two  extreme 
ends  of  the  panel  system  are  supported  by  stringers.  In  this  case  we 
write 


HD 


where  the  state  vectors  are  for  the  free  edges  just  outside  the  lines  of 
attachment  of  the  end  stringers.  Hence  in  both  state  vectors  the 
moment  and  shear  are  zero;  ie. 


i 


0 


'N 


0 


■SO 


Thus  the  following  two  by  one  matrix  equation  can  be  extracted  from 
equation  (41). 


r 

—  — 

t31  t32 

l 

l 

no 

= 

* 

> 

N 

Hi  tU2_ 

0 

Y  ' 

^  n4> 

(1*2) 


For  this  equation  to  have  a  non-trivial  solution,  we  must  have  the  fre¬ 
quency  determinant  equation 


N 


31 

t32 

41 

%k2 

=  i(w)  ®  0 


(1.3) 
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which  can  be  solved  for  the  natural  frequencies  of  the  skin-6tringer 
system.  The  simplist  method  of  obtaining  a  numerical  solution  to 
this  complicated  transcendental  equation  is  that  of  graphing  A(co).  Its 
zeros  are,  of  course,  the  natural  frequencies. 

For  a  given  root  to  of  the  frequency  determinant  equation,  the 
associated  normal  mode  can  be  obtained  as  follows.  Let  ^at  ( l ,  0) 
be  arbitrary.  Then  the  magnitude  of  at(^*0)  is  determined  from 
equation  (42)  as 


n 


Let  us  say  we  are  interested  in  the  modal  deflections  of  the  j^h  panel. 
By  use  of  transfer  matrices  we  can  obtain 


(WO 


(*»5) 


The  elements  of  { Y } j are  the  constants  a0  »  Ai»A2*  andA3 


appearing  in  equation  (27a)  from  which  we  can  calculate  Yn(y). 
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VI.  CURVED  STRINGER -PANEL  SYSTEMS 


The  analysis  for  curved  panels  is  more  involved  that  that  for 
flat  panels.  It  is  possible  to  approach  the  problem  by  means  of  a 
distributed  mass  mathematical  model,  but  to  reduce  the  complexity  of 
the  problem  we  will  lump  the  distributed  panel  mass  along  discrete 
"mass  lines'*  running  parallel  to  the  supporting  stringers.  The  mass 
lines  will  be  numbered  .  .  .  ,  j-1,  j,  j+1,  ...  .  The  corresponding  masses 
per  unit  length  of  mass  line  are  designated  ...»  5  >  ^j+1  ’  ‘ ’ 

and  the  rotational  moments  of  inertia  per  unit  length  of  mass  line  are 
•••j  Jj-1’  *^j+l5  •••  •  The  line  masses  are  connected  by  massless 

segments  of  thin  cylindrical  shells.  This  discrete  approximation  of  a 
curved  panel  is  depicted  in  Figure  11. 

The  usefulness  of  this  model  is  as  before  dependent  upon  being 
able  to  consider  the  panels  as  simply  supported  at  the  frames,  i.e.  along 
the  curved  edges.  When  this  ability  exists  the  problem  is  reducible  to 
one  of  one  dimension,  and  the  same  general  procedure  as  carried  out  in 
our  previous  problems  will  yield  a  solution. 


We  will  first  concern  ourselves  with  developing  the  field  transfer 
matrix,  the  matrix  which  will  carry  us  across  the  massless  shell  seg¬ 
ments  between  mass  lines.  We  will  base  our  development  of  the  field 
transfer  matrix  on  Donnell1  s  shell  theory,  an  approximate  theory  which 
is  obtained  by  eliminating  those  terms  of  the  more  exact  theory  which  are 
of  minor  importance  to  the  case  at  hand.  For  an  unloaded  and  massless 
cylindrical  shell  element,  the  Donnell  equations  tU  are 


2  32u 
a  - 

a  2 
dx 


+  A  (1-vJ  +  —  a(l+v) 

2  a*2  2 


32\ 


3x3<J> 


■va 


3w 

3x 


0 


j  a  (1+v) 


32 


u 


3x3<{> 


1 

2 


a2(i_v)i^ 

3x2 


=  0 
3<J> 


(U6) 


va 


3u 

3x 


+ 


3 v 
3<f> 


w 


&4  94v  +2a2  d4w 

3x4  3X2  Y  4*2 


0 


where 
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u 


Figure  11.  A  Typical  Curved  Panel 
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the  radius  of  curvature  of  the  segment  (note  that 
this  quantity  may  change  from  segment  to  segment) 


k  =  h2/12a2 

v  =  Poisson's  ratio 

<p  =  angular  coordinate  of  the  shell  segment 

(a<P  corresponds  to  y  of  the  flat  panel) 

These  are  two  second  order  and  one  fourth  order  partial  differential 
equations  in  three  unknowns,  the  displacements  u,  v,  and  w.  They  can 
be  combined  into  the  following  single  equation  of  the  eighth  order  in 
one  unknown. 


a8  kV8  w  +  a4(l-v2)  —  —  =  0 

3x4 


(UT) 


and  V8  =  (?2)4 


With  simple  support  conditions  at  the  frames,  the  solutions  for  the  case 
of  a  single  frequency  oscillation  can  be  written  as 


u 


v 


(M) 


w 


e 


iwt  $n(4,)  sin£I2L 
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Substituting  the  above  expression  for  w  into  equation  (47),  we  find  that 
#n(4')  must  satisfy  the  following  ordinary  differential  equation. 


(8) 


-  k 


2 

Si 


*(6) 

n 


6q£  *(U) 


n 


* 


i  t 

n 


where 


Si  = 


+  <  K  +  -  0 

V12a2(l-v 
h2 


(49) 


and  p 


-v2) 


The  eight  characteristic  roots  of  this  equation  can  be  written  as 


or,  to  better  suit  our  purposes,  in  real  and  imaginary  component  form 
a  s 

-  (y1  -  i^)  and 

±  (y  2  1 


The  values  of  these  components  work  out  to  be 

Y  =  /T"( q.2  +  J2  a  p  +  p2)1^  cos(^-  arc  tan  - P -  ) 

1  «  n  ^  *  ^qn  +  P 

61  -  +  &  v  +  p2)lA  sin(l arc  tan  ] 

Y2  -  J2  qnp  +  p2)1^  cos[|-  arc  tan  p/(/2  qn  - 

62  ~  ~  S2  %P  +  P2)1^  sin[|-  arc  tan  p/('/2  qR  - 

and,  in  all  cases  the  principal  value  of  the  arctangents  are  to  be  used. 

When  the  characteristic  roots  are  expressed  in  terms  of  their  real  and 
imaginary  components,  we  can  say  that  Eq.  (49)  is  satisfied  by  any  of 
the  following  eight  beharmonic  functions 


(50) 

P)] 

P)] 
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coshY-^4>  * 

cosS^<t> , 

sinhy^(j)  * 

sin 

coslry^  * 

cos624>. 

sinhy2(t)  * 

sin 

624>, 

sinhy^  • 

cos6^<j>. 

coshy^  # 

sin 

sinhy^^  * 

cos624). 

coshY24)  # 

sin 

*5  2  4*  • 

The  eight  solutions  we  will  use  will  be  linear  combinations  of  the  above 
solutions  so  as  to  facilitate  the  construction  of  the  field  transfer  matrix. 
We  will  write  our  desired  solution  in  the  form 

7 

*(♦)  =  £  A.  f.U)  (51) 

n  i=0  1  1 

and,  so  as  to  have  even  subscripts  of  f  imply  even  functions  of  <f> ,  we 
further  write  for  i=0,  2,  4,  and  6 

f\  =  a^coshY^  cos 6^  +  ou0  sinhY^<j>  sin  6^ 4> 

(52) 

+  ai3coshY2(tl  00862,1,  +  sinh^  sin 

and  to  have  odd  subscripts  of  f  imply  odd  functions  of  <j> ,  for  i=l,  3,  5  and 
7 


a^sinh  cos  6-^ 

+  a^2  cosh 

Y1<(> 

sin 

oii3  sinhY^  cos 

+  cosh 

Y24> 

sin 

62<)> 

(53) 


Before  we  can  construct  a  field  transfer  matrix,  we  must  agree  on 
a  suitable  state  vector.  The  eight  quantities  necessary  to  specify  the  state 
of  the  cylindrical  shell  segment  are 


7*'  ,  (M.)  ,  (V  )  (N.)  .  and  (N  )  . 

a  n  <j>  n  4>n  <|>  n  <?>x  n 


These  eight  quantities  in  the  order  given  will  constitute  our  state 
vector.  The  meaning  of  the  last  four  quantities  is  to  be  seen  from 
the  following  relations t ? 1 .  The  second  part  of  these  relations  are 
part  of  Donnell* s  theory  of  cylindrical  shells. 


M,* 


=  e] 


(M4>)n 
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-  .l“t  (V.)n  sin  up  -  D  [  4  .  +  O-v)  -1 


3  4)' 


33w 


e1(i,t  00  n  sin  ^2L  =  -Jh^  i  iv  .  w  +  y  V \ 
z  i-v2  \  a  3<j»  a  3X  / 


and 


>x 


=  e 


iwt 


/lT  \  -  T1TFX 

^♦x^n  sin  o 


Eh 


2 (1+v) 


/I  . 

y  a  3<J>  3x  J 


(54) 


See  Figure  12  for  visualization  of  M$,  V<f>,  N$,  and  N$x.  From  these 
relations,  Eqs.  (54),  and  their  counterparts  for  the  deflections, 
Eqs .  (48) ,  we  can  see  that 


<Vn 

depends  upon 

a" 

n 

and 

*n 

<Vn 

depends  upon 

n 

and 

*'« 

<yn 

depends  upon 

*  n* 

V 

and 

<V>» 

depends  upon 

T  * 

1  n 

and 

Therefore  our  initially  chosen  state  vector  can  be  replaced  by  a  state  vector 
composed  of  the  eight  quantities  <t>n»  $'n,  'n» 


v '  T  and 
T  n  n 


T'  Moreover,  from  the  second  and  third  expressions  of  Equations 


we  find  that  Tn  can  be  expressed  in  terms  of  f'n 


T'n,  in  terms  of  $'n 


and  4” 


Thus 


(in  4'n 


*  If 
n» 


and  *n 


*"n 


*n(4> 


*n  rn  !^J  is  still  another  possible  state  vector.  It  is  this  latter 

form  of  the  state  vector  that  we  shall  use  as  a  convenient  temporary 
replacement  for  our  initially  chosen,  more  descriptive  state  vector  in 
order  to  determine  the  essential  features  of  the  field  transfer  matrix. 
Once  this  is  done  we  will  use  a  conversion  matrix  to  return  to  our 


initial  state  vector. 


At  this  point  we  return  to  equation  (51)  and  resume  our 
search  for  the  field  transfer  matrix.  Another  attribute  that  we  will 
require  of  our  solution  is  that  the  f^($)  are  chosen  so  that 
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their  Positive  Directions 


A0  =  $n(0) 

> 

H 

It 

1©. 

D  ~ 

O 

a2  =  *"n(0) 

a3  =  *™(0) 

A*  =  ^(4)W 

S  ■  »n(°) 

(55) 

O 

EM 

II 

C 

A?  =  T'^CO) 

One  further  consideration  before  we  fully  develop  our  solution:  the 
solution  for  $  n  is  not  independent  of  *n.  From  the  second  and  third 
equations  of  (46)  we  may  obtain 


r  + 
n 


'“In2  *n 


Vn  + 


n2<t> 


in  + 

n 


(5) 


'3  n 


(56) 


where 


l-v 


=  -2  (  )  k  qn2 


l-v 


n3 


1+y 

l-v 


Now  let  us  define 

7 

¥n(4>)  =  I  A i  SiU) 

i=0 

Then 

7  7  ( 5) . 

I  tiAe\  +  vq^g.)  =  1  Ai(nx  f'i  +  n2  f'"i  +  n3  V  > 

1  1  n  1  i=0 


i=0 


(57) 


Since  this  is  true  for  any  set  of  A  - *  Eq.  (56)  becomes  a  relation  be¬ 
tween  the  functions  g±( <J>)  and  ^(4) for  any  particular  i;  i.e. 


vq  2gi  =  fV  +  npf'i  +  TI3 


(5) 


(58) 


We  note  that  { <(> )  is  an  odd  function  if  fj(<f>)  is  an  even  function,  and 
vice  versa.  A  property  of  the  f  j  functions,  which  we  will  use  again 
later,  is  their  essential  permanence  of  form  under  differentiation  ex- 
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cept  that  the  odd  functions  become  even,  and  vice  versa.  Thus  a 
for  g^(<J>)  may  be  written  in  the  form 


6iU) 

— 

8il 

sinh 

Yl4> 

COS 

6±<t, 

+ 

B12  cosh  yi4> 

sin  514> 

+ 

6i3 

sinh 

Y2 ♦ 

COS 

62<t> 

+ 

8il*  cosh  Y2<t> 

sin  62<p 

for 

i=0 

,2,4 

,  6;and 

g.U) 

= 

8il 

cosh 

Yi4> 

cos 

6^4) 

+ 

6i2  sinh  <j> 

sin  .  6^4> 

+ 

6i3 

cosh 

Y2<t> 

cos 

62<p 

+ 

eilt  sinh  y2<p 

sin  62$ 

for  i=l,  3,  5,  and  7. 

It  will  be  convenient  to  collect  these  various  definitions  of  ^and 
into  matrix  form  as: 
for  i=0,  2,  4,  or  6. 


f^U)  = 

bu 

ai2 

ai3  1 

“uj 

Ue) 

ei(4>)  = 

L6ii 

&i2 

6i3 

Bi*J 

Uo> 

for  i=l,  3, 

5,  01 

•  7 

f±(*)  = 

|_?il 

ai2 

“13  “Ikj 

U0} 

ei(<l>)  = 

L%i 

$i2 

6i3  bH.J 

ue} 

where 

Ue>  = 

f 

cosh 

Y2<t> 

cos  6-^4>  ^ 

sinh 

Y24> 

sin  6^4> 

r 

cosh 

Y2<<> 

cos  6^4) 

sinh 

Y2<t> 

sin  6^4) 

J 

solution 


( 59a) 


8i 


( 60a) 


(60b) 
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{Cq} 


sinh 

cos 

cosh 

Yl<f> 

sin 

Sl<|) 

sinh 

Y2 ♦ 

cos 

cosh 

72$ 

sin 

J 


To  obtain  our  desired  solutions  for  and  we  have  to  evaluate  four 
constants  for  each  of  eight  f^  functions  and  each  of  eight  functions 
or  a  total  of  sixty-four  constants.  We  begin  by  returning  to  our  differ¬ 
ential  relation  between  $>n  and  4^  in  order  to  relate  the  6  coefficients 
to  the  a  coefficients,  and  thereby  reduce  immediately  the  problem  to  that 
of  evaluating  thirty-two  coefficients. 


As  a  step  preliminary  to  shaping  the  relation  between  the  a  and  8 
coefficients,  we  want  an  economical  means  of  describing  differentiation 
When  properly  organized,  straight  forward  differentiation  yields 


=  L?il  °i2  ai3  °iiJCD] 

for  l=even ,  and 


_d_ 

d* 


fiU)  =  L°il  ai2  ai3  °ibJ  ^ 


for  i=  odd,  where 


ID1 


-6n 


Yi 

0 

0 


(6la) 


(61b) 


Of  course,  the  similar  situation  that  exists  for  gi (<|>)  can  easily  be  deduced 
from  the  above.  Therefore  differentiation  of  either  or  g^)  is 

equivalent  to  the  application  of  the  matrix  of  transformation  [D]  to  the 
opposite  column  matrix  of  biharmonic  functions;  that  is,  for  example, 
differentiation  of  with  i  even,  means  replacing  {£g}  with[D]UQ}. 


Now  we  are  ready  to  relate  the  8  coefficients  to  the  a  coefficients 
by  means  of  Eq.  (58).  Direct  application  of  the  differential  operator 
yields 


LBiJ  (  ID]  +  vx^2  [I]) 


=  J (^[D]  +  n2  [Dl  +  n3LDls) 


49 


or 

U.JiK]  *  LOjJIl] 


UiJ  = 


Oi[L][K] 


-1 


(62) 


where  (I]  is  the  identity  matrix,  and  the  definitions  of  fK]  and  f L 1 
are  clear.  Let  the  coefficients  for  g'<$)  be  k„,  Ki9,  kh  ,  and  , 
and  those  for  W  be  ,a,  ri2,  T^/Hhen2’  i3’  i4’ 

UiJ  =  Uj  CD]  [L]  [K]  "*1 


biJ  =  LaJ  CD]  [L]  [K] 


-1 


The  task  that  remains  is  that  of  specifying  the  thirty-two  a 
coefficients.  We  do  this  by  returning  our  attention  to  Eq.  (5l), 
Eqs.(52)  and  (53),  and  Eq.  (55)-  It  is  our  purpose  to  guarantee 
that  Eq.  (55)  satisfied.  Since  (l)  we  segregated  the  functions 
f£($)  and  g^(t}>)  into  odd  and  even  functions  ,  since  (2)  we  are  inter¬ 
ested  in  the  values  of  these  functions  for  zero  arguments,  since  (3) 
and  even  number  of  derivatives  of  odd  functions  and  an  odd  number  of 
derivatives  of  even  functions  produce  odd  functions,  and  since  (^) 
the  value  of  any  odd  function  of  zero  argument  is  zero,  we  may  con¬ 
veniently  treat  the  odd  and  even  functions  separately  as  follows. 

The  sixteen  conditions  the  even  functions  must  satisfy  in  order  to 


fulfill  Eqs ,  (55) 

are  concisely 

stated  in  the 

f0(o) 

f0"(0) 

f0(l4)(o) 

s'o<°>- 

f2(0) 

f2”(0) 

f2(4)(o) 

8'2<0) 

*1,(0) 

V(0) 

*4^(0) 

s'l,(0) 

*6(0) 

f6"(0) 

f6(lt)(0) 

s'6(o) 

[I] 


(63a) 


The  sixteen  conditions  for  the  odd  functions  are 
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[I] 


<  63~b) 


~t\(o) 

gi(0) 

s\(0 T 

f’3(0) 

f"'3(0) 

g3(°) 

g"3(0) 

f’5(0) 

f»«5(0) 

g5(°) 

g"5(0) 

f"’7(0) 

gT(0) 

g"7(0) 

One  advantage  of  treating  the  even  and  odd  functions  separately  is  that 
it  reduces  the  size  of  the  above  matrix  statements.  This  is  a  total  of 
thirty- two  conditions  for  the  thirty-two  a  coefficients.  Let  us  look  at 
how  these  statements  in  terms  of  the  f ±  (+)  and  gj*)  functions  may  be 
conveniently  rephrased  in  terms  of  the  a  coefficients.  Again  considering 
the  even  functions  as  examples,  fo(0>  =  1  means 


f0"(0)  =  0  means 


LoqJ  [3)]^ 


=  l 


=  0 


f0(“)(0)  =  0  means 


i 

.0. 


La0J  [»)**  < 

g’0(0)  =  0  means 

L“qJ  [D][L][K]_1  J  0 


0 


0 


If  we  repeat  this  analysis  for  the  remaining  twelve  expressions,  and 
combine  our  results,  they  are 


“oi 

a02 

“03 

e 

o 

-P-. 

a21 

“22 

“23 

“24 

°Ui 

“42 

“43 

“44 

“61 

“62 

“63 

“64 

(64) 
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where  [Q]  has  as  its  columns 


1 — 1 

V 

f  \ 
1 

0 

[»]2^ 

0 

► 

i 

1 

0 

J 

0 

J 

'i 

1 


if’ 


Hence,  for  the  even  functions 


and  t»I[L][K]_1  « 


0 


[A]  =  [QpX 


(64a) 


and  for  the  odd  functions 


[A]  =  [Up1 

where  [U]  has  as  its  columns 


m  i 


(64b) 


r  y 
I 

ri 

1 

1 

V 

0 

0 

0 

0 

A 

[Dl3  < 

>  ,  [L]  [K]_1  -< 

and  [3)]2[LJ[K]-1  < 

1 

1 

1 

1  ' 

0 

0 

0 

o  ; 

v  y 

V.  J 

J 

At  this  point  we  have  determined  all  the  functions  p($>)  and  g.($)  in 


*n(<t>) 


I  At  ti  <« 

i=0  L 


(51) 


V« 


I  A  g  («j>) 

i=0 


(57) 


and  we  did  this  so  that  Eqs.  (55)  are  satisfied.  Hence  the  state  vector 

,(«) 


LsJ  =  L‘n  *'n  *i'  *"  *»  ’  *a  *i  *i’J 


(65) 
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can  be  determined  at  any  value  of  the  argument  <(>  from  its  value  at 
=  o.  That  is,  in  particular  we  now  have  the  transfer  matrix  by 
which  to  obtain  {5}  at  <t>  =  0^  . 


{  =  } 


£ 

J 


{5} 


r 

J-l 


where  the  elements  of  [r]j  are  given  by 

rk*=|_ail  °i2  °i3  a14j[D]  k-l{e(e  )} 

for  l  =  itl,  k  =  1,  2,  3,  4,  and  5 


and 


rk£  =  Lail  ai2  °i3  «iU  J  [n3k“6  [L]  [K]~±  ( C  (6  j)  ) 


-1 


for  £  =  i+1 ,  k  =  6 ,  7 ,  and  8 

and  where 

U(0j)>  =  Ue(0j)}  for  £+k  =  even 

U(0,)}  =  U  (0,))  for  £+k  =  odd 

J  o  J 

The  {E}  vector  is  not  suitable  for  the  incorporation  of  boundry  con¬ 
ditions,  nor  does  it  conveniently  mesh  with  our  analysis  of  stringers. 
Therefore  we  will  now  convert  LHjto  LZj  =  [T n  H'n  4>n  C )n ( >n ( >n ( ) 

The  first  component  of  our  mixed  displacement  and  force  type 
state  vector  can  be  expressed  in  terms  of  the  components  of  L-  J  ^y 
referring  to  the  third  of  Eqs.  (46).  From  this  we  find  that 


Tn($)  =  —  [r  - 

11  v  n 


(l+kq£) 


<t>  +  2kq2  <t> 


n 


The  second,  third,  and  fourth  elements  of  LZJ  present  no  difficulty. 

Referring  to  Eqs.  (5*0,  we  have 


M<j,  = 


liiL  +  ...  g2w  ) 

3  <t>2  ix2  1 
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So, 


<Vn  *  2^*"n  ♦> 

az 


Similarly 


-  a  > 


and 


(Vn  =  T[*»'  -  (2-',)qnt'n' 

a3 

Similarly 

N  =  I  i  lx.  -  2.  +  V  iH.\ 

<J>  i_v2  V  a  d<f>  a  3X  y 


so 


(Nj 


Eh 


n  a(l-v2) 


Similarly 


-*n-v%  Tn> 
m  ^  ^qn*n  “  2qn*"n  + 


„  _  Eh  ,  1  3u  .  3v  x 

*x  2(l+v)  (  a  34>  +  9x 


SO 


(N*x)n  “  2a(l+v )  (T'n  +  V 

It  follows  from  the  second  of  Donnell's  equations,  Eqs.  (46)  that 

T’n  =  (Hv  )qn  *"n  "  lT7  qn  "  (i+v)q  *n 


Therefore 


(N.  )  =  — — (  —  T"  +  vq  V  -  *'  ) 


*x'  "  a(l+v)2  q 


n 


n 
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Now  we  can  set  up  our  conversion  matrix  so  that 


{z}  =  [B]j  {5}  (66) 

where  the  subscript  j  is  associated  with  [b]  indicates  that  this  matrix 
is  to  be  computed  from  the  physical  constants  of  the  j**1  shell  segment. 
The  details  of  the  matrix  [b)  j  can  be  found  on  the  following  two  pages. 

The  {Z}  vectors  at  the  two  opposite  edges  of  the  jth  shell  segment  can 
now  be  related  by  writing 


U>;  -  [Blj  I*la  [B]-1 


(6T) 


This  completes  our  development  of  the  necessary  field  transfer 
matrix,  and  we  now  turn  our  attention  to  the  point  transfer  matrix 
mecessary  to  account  for  the  effects  produced  by  a  stringer  and  a 
mass  line.  Eqs.  (37)  are  our  starting  point  for  the  case  of  a  stringer. 
If  we  insert  the  series  expansions  for  the  deflections,  monent,  shear, 
and  tensile  force  into  these  equations,  we  can  write  the  result  for  a 
stringer  as 


(M  )r  -  (M.)* 
9  n  9  n 

<Vn  -  'Vn 


c  -  $’  +  a  $  +  f  <r 

a  n  n 


n 


(68) 


(H.)r  -  (N.  P  =  h  T  + 
'  9  n  9  v  n 


£  i 


+  g 

n  n 


where  the  quantities  c  ,  d,  and  e  are  detailed  in  Eqs.  (39),  and 
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f  =  EIc  sz  (nr>4f  pa(cz  -  sz>“2 

g  =  EI  /M.\*  (69) 

nc''  £  ' 

h  =  EI?  (^L)i  pAu>2 


The  corresponding  contributions  of  a  panel  skin  mass  line  located  at 
the  same  station  are-tTw2 ,  +pw2 ,  and  -pmj2  respectively.  Thus  we  may 
compose  a  general  point  transfer  matrix  as 


1 


0  0 
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0 
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Of  course,  if  no  stringer  exists  at  station  j,  we  merely  set 

A  A  A  A  A  A 

c  =  d=:e  =  f=  g  =  h  =  G,  On  the  other  hand*  if  there  is  no  panel 
mass  line,  pj  =  Jj  =0* 

Now  that  we  have  described  the  basic  transfer  matrices  [F]  and  [G], 
we  can  relate  the  state  vectors  at  any  two  locations  in  the  panel 
system.  We  determine  the  natural  frequencies  as  before  by  relating 
the  state  vectors  at  the  two  ends  of  the  panel  row.  For  example,  if 
both  ends  of  the  panel  row  are  simply  supported  (i.e,  =  =  (M^)^ 

3  (N(j)n  =  0),  these  two  vectors  and  the  connection  between  them  are, 
for  any  value  of  n 
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0 

v 

0 
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l  r 

where  [T] 

N  0 
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=  [F]n  [G] 
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SL  r 
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‘t’X 
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[G]1  [F]r 


Then  proceding  as  before,  the  determinant  equation  for  the  natural  fre¬ 
quencies  is,  for  any  n, 
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C24 

C26 

t28 

31 

C36 

C38 

51 

C54 

C56 

C58 

71 

C74 

C76 

C78 
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=  A(u)  =  0 


0 


(70) 


The  determination  of  the  normal  mode  associated  with  a  natural  frequency 
can  in  theory  be  approached  in  a  manner  analogous  to  the  case  of  flat 
panels.  However,  in  both  the  flat  and  curved  panel  cases,  difficulties 
may  arise  when  actually  carrying  out  these  frequency  and  mode  shape 
computations . 
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VII.  NUMERICAL  COMPUTATIONS 


Two  types  of  difficulties  can  arise  when  carrying  out  the  numerical 
solution  of  any  of  the  frequency  determinant  equations.  Reference  [2], 
section  7.1  is  a  good  and  sufficient  explanation  of  these  difficulties. 

The  first  of  these  difficulties  is  the  loss  of  significant  figures  when  the 
frequency  determinant  is  reduced  in  an  ordinary  manner.  This  case  is 
exemplified  by  Eq.  7-3  of  Reference  [2].  The  second  difficulty  arises 
when  the  elastic  supports  are  very  stiff  in  comparison  with  the  remainder 
of  the  structure.  The  error  introduced  is  explained  by  Eq.  7-^.  In  sec¬ 
tions  7-2  and  7-3  there  are  presented  two  methods  of  avoiding  these 
obstacles.  The  first  method,  the  Delta  Matrix  method,  is  limited  in  the 
size  of  the  frequency  determinant  to  which  it  can  be  applied.  It  is  not 
difficult  to  apply  successfully  to  a  four  by  four  frequency  determinant, 
such  as  that  of  a  flat  panel  system,  and  by  dint  of  much  programming, 
it  can  also  be  applied  to  an  eight  by  eight  frequency  determinant .  A 
larger  determinant  would  probably  be  beyond  the  capabilities  of  most 
present  day  digital  computors.  In  section  7-^  a  trial  and  error 
scheme  is  explained.  A  short  summary  of  this  method  is  available  in 
Reference  [l],  section  IV. 

Even  when  a  direct  graphical  solution  of  the  frequency  determi¬ 
nant  equation  is  no  problem,  there  is  a  possibility  that  the  values  of 
A(u))  may  pass  through  zero  very  abruptly!*?].  Thus,  while  it  is  easy  to 
make  a  very  good  estimate  of  w  per  se,  it  is  very  difficult  to  obtain 
near  zero  values  of  A(a))  for  tne  purpose  of  computing  the  mode  shapes. 

As  a  result,  if  many  transfer  matrices  are  necessary  to  describe  the 
structure,  the  error  in  the  initial  state  vector  will  grow  to  the  point 
where  the  final  state  vector  is  not  even  approximately  representative  of 
the  given  boundry  condition  at  that  end.  A  means  of  circumventing  this 
error  accumulation  is  that  of  adopting  a  forced  vibration  point  of  view. 

To  begin,  sinusoidal  shear  forces  (or  moments)  of  frequencies  just  on 
either  side  of  the  natural  frequency  of  interest  are  applied  in  turn 
to  the  structure,  and  the  respective  states  of  the  structure  are  calcu¬ 
lated.  If  it  is  found,  as  would  be  expected,  that  the  computed  states 
corresponding  to  the  two  such  frequencies  are  very  nearly  the  same,  then  it 
may  be  concluded  that  both  states  are  dominated  by  the  same  mode — the  mode 
corresponding  to  the  bracketed  natural  frequency.  Therefore  it  is  reason¬ 
able  to  call  the  averages  of  the  deflection  components  of  these  states  the 
mode  shapes  of  that  nearby  natural  frequency. 

The  following  example  will  illustrate  the  features  of  the  forced 
vibration  procedure  of  calculating  modal  states.  Let  us  apply  a  sinusoidal 
shear  loading  of  a  unit  amplitude  at  station  s  of  a  curved  panel  with 
simple  support  end  conditions.  By  combining  our  successive  transfer 
matrix  equations  we  arrive  at 
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[T]J  <Z>0  +  i[< 


"\ 


where  the  transfer  matrices  are  evaluated  at  one  of  the  two  selected 
frequencies.  Again  we  make  use  of  the  boundry  conditions  and  reduce 
the  above  eight  by  one  matrix  equation  to  the  following  four  by  one 
equation 
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This  matrix  equation  can  be  solved  for  LT  *'/a  V  N<j,  Jq 
by  inverting  the  four  by  four  reduced  transfer  matrix.  This  inversion 
is  always  possible  even  in  the  case  of  an  undamped  model  analysis  be¬ 
cause  the  determinant  of  this  matrix  is  identical  to  A(w)  which,  of 
course,  is  non- zero  by  our  choice  of  a  frequency  other  than  a  natural 
frequency.  Furthermore,  the  more  abrupt  the  change  in  the  value  of 
A(w)>  ije.  the  more  acute  the  difficulty,  the  closer  the  values  of  the  selected 
frequencies  can  be  to  the  natural  frequency.  (Differences  of  0.  1  radians 
have  been  satisfactory  in  all  respects  in  the  cases  known  to  the  author.) 

It  is  of  course  necessary  to  avoid  nodal  points  when  applying  the  driving 
force  or  moment.  This  can  be  assured  by  altering  the  point  of  application 
of  the  driving  load  and  then  comparing  the  results. 
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followed  by  an  explanation  of  the  application  of  the  transfer  matrix 
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